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Abstract 

The Podolsky generalized electrodynamics with higher derivatives 
is formulated in the first-order formalism. The first-order relativistic 
wave equation in the 20-dimensional matrix form is derived. We prove 
that the matrices of the equation obey the Petiau-Duffin-Kemmer al- 
gebra. The Hermitianizing matrix and Lagrangian in the first-order 
formalism are given. The projection operators extracting solutions of 
field equations for states with definite energy-momentum and spin pro- 
jections are obtained, and we find the density matrix for the massive 
state. The 13 x 13-matrix Schrodinger form of the equation is derived, 
and the Hamiltonian is obtained. Projection operators extracting the 
physical eigenvalues of the Hamiltonian are found. 

1 Introduction 

There is currently a renewal of interest in higher derivative (HD) field the- 
ories. HD field equations appear in many models such as renormalizable 
quantum gravity [I], Podolsky's generalized electrodynamics j2], the Lee- 
Wick model [2] and others. One of the reasons to consider HD theories is 
to improve renormalization properties of theories and to remove ultraviolet 
divergences [I]. However, HD models suffer some difficulties connected with 
the presence of ghosts. These can lead to the violation of unitarity [5], [6]. 
Nevertheless, in some HD models these problems with negative probabilities 
and S-matrix unitarity can be avoided [7]. Also, the quadratic divergence 
associated with the Higgs mass were removed in the HD Lee- Wick standard 
model [8], that solves the hierarchy problem. Extensions of the minimal 
standard model living to new physics are justified until observations at the 
Large Hadron Collider (LHC) will be analyzed. 

It is well known that in classical electrodynamics the electromagnetic 
mass is infinite and, therefore, there are infinities associated with a point 
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particle. One of the ways to solve this problem in classical theory is to 
use the Lorentz invariant regularization of the Maxwell equations at short 
distances. With the help of an appropriate cutoff the point particle limit can 
be achieved. This programme was realized in Podolsky's electrodynamics. 
Firstly the interest to Podolsky's electrodynamics was due to the fmiteness 
of the theory: the electromagnetic energy of a point charge is finite contrarily 
to ordinary electrodynamics. If distances are much greater than a cutoff, 
Podolsky's electrodynamics converts into Maxwell's electrodynamics. The 
solution to Poisson equation for the potential corresponding to a point charge 
e, located at the origin, in Podolsky's electrostatics is given by [2] 



where a is a new parameter of the theory with the dimension of the length 
playing the role of the cutoff. This potential becomes the Coulomb potential 
at distances much bigger than a. At r — > the potential <p approaches the 
finite value e/47ra. The energy of the field for a point charge is also finite 
in the hole space. Thus, the electrostatic energy can be considered as the 
regularized electromagnetic mass of a point charge. It was shown [9] that 
higher derivatives terms in Podolsky's equations suppress unphysical runaway 
solutions with exponentially growing acceleration of the Abraham-Lorentz 
equation. There are not unwanted solutions if the cutoff is greater than half 
of the electron classical radius. The upper bound on the parameter is of the 
order a ~ 10~ 16 cm [9], i.e. the same as the Compton wavelength of the 
neutral Z-boson. Classical Maxwell's electrodynamics is not valid at small 
distances and time intervals due to quantum effects. It was also mentioned 
in [TU] that in the framework of non-relativistic quantum theory a natural 
cutoff of order of the electron Compton wavelength is effectively appeared by 
QED processes in close analogy with the classical theory of extended charges. 
Thus, one may treat the classical Podolsky's electrodynamics as an effective 
theory where a cutoff introduced, a, is due to the quantum processes at small 
distances (large momentum). If distances are larger than a, the classical 
regime begins. 

At the same time although QED describes all experimental data well, 
there are some internal difficulties with the regularization [11]. We mention 
infrared catastrophe: when the average number of photons n — > oo, then the 
matrix element | < out\0 in > | — >• 0, and it is impossible to construct 
the "out" Fock space from the "in" space, nor to find unitary operator S 
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[TT] . The authors [12] wrote: " There is an alternative possibility to avoid in- 
frared divergences. We give the photon a small mass \i. This will cut off the 
low- energy region since now k° > \x and therefore remove the infrared diver- 
gence. . . . The infrared divergences of quantum electrodynamics are essentially 
classical" . We continue with the citation [13] : " Another aspect of infrared 
singularities related to the long-range character of the Coulomb forces. The 
latter induces an infinite phase shift on the scattered plane waves. To prevent 
it we may introduce a screening factor which in a consistent theory would 
be related to the fictitious photon mass //" . Thus, in QED the cutoff is intro- 
duced "by hands" as for small distances (to remove ultraviolet divergences) 
as well as for large distances (to avoid infrared catastrophe). Therefore, one 
may consider naturally to extend classical Podolsky's electrodynamics on 
the quantum level where the cutoff is appeared due to the presence of higher 
derivatives. Anyway, different aspects of Podolsky's electrodynamics, in our 
opinion, have a definite theoretical interest. 

Some features of Podolsky's electrodynamics were investigated in [13] , 
[To] . [TB] . [T7] . The goal of this paper is to formulate Podolsky's electrody- 
namics equation in the form of the first-order relativistic wave equation, and 
to obtain solutions in the form of projection matrices. 

The paper is organized as follows. In Sec. 2, the third-order field equa- 
tion is discussed. We derive the first-order relativistic wave equation for 
Podolsky's electrodynamics in the 20-dimensional matrix form. The Hermi- 
tianizing matrix and the Lagrangian in the matrix form are found in Sec. 3. 
The projection operators extracting solutions of field equations for definite 
energy and spin states of particles are obtained in Sec. 4. We find the density 
matrix for the massive state. In Sec. 5 the 13 x 13-matrix Schrodinger form 
of the equation is derived, and the Hamiltonian is obtained. Solutions of 
this equation are found in the form of projection operators. The results are 
discussed in Sec. 6. In Appendix A, we consider the first-order wave equa- 
tion in the presence of the charge current density. The Lorentz covariance 
of the equation is proven. Some useful products of matrices are derived in 
Appendix B. We obtain "minimal" polynomials of the matrix of the equation 
for massless and massive states. In Appendix C the "minimal" polynomial 
of the Hamiltonian matrix is derived. 

The Heaviside's units are chosen, h — c — 1, and Euclidian metric is used, 
Xfi = {xm, ixo). Greek letters range from 1 to 4 and Latin letters range from 
1 to 3, and there is a summation on repeated indexes. 
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2 Field equations 



2.1 Third-order field equations 

The Lagrangian of Podolsky's electrodynamics is given by [2] 



lp2 ,.2 i ,2 

2 



(1) 



where = d u A^ is the field strength, d u = d/dx u = (d/dx m , d/d(it)). 

The dimensional parameter a can be written as a = 1/m, where m is the 
mass parameter. The Euler-Lagrange equations follow from Eq.(l): 



(d 2 a - m 2 ) d,F u , = 0. (2) 

The Lagrangian (1) and the equation of motion (1) are gage- invariant under 
the [/(l)-group. We can represent Eq.(2), in the momentum space as the 
matrix equation: 

(p 2 + m 2 ) (p 2 - p ■ p) A = 0, (3) 

were A = {^4^}, the matrix-dyad p-p, with matrix elements {p-p)^ v = V^Pvi 
is introduced, and the four- momentum being p^ = (p, ipo)- The matrix 
M = p 2 — p ■ p obeys the minimal polynomial M(M — p 2 ) = 0, so that 
the eigenvalues of the matrix M are zero and p 2 . Thus, Eq.(3) leads to the 
dispersion equation: 

p 2 (p 2 + m 2 ) = 0. (4) 

Eq.(4) shows that there are massless and massive states in the spectrum. 
The propagator of fields is given by 

1 1 (5) 



p2 (jp. _|_ m 2^ p2 p2 _|_ m 2 

The first term in Eq.(5) is the propagator of the photon massless field and 
the second term corresponds to the propagator of the massive state of the 
field. A "wrong" sign (— ) in Eq.(5) indicates that the massive field state is 
a ghost. As a result, the massive field state gives the negative contribution 
to the energy [2J, and the classical Hamiltonian is unbounded. To have 
the positive eigenvalues of the Hamiltonian in the second quantized theory, 
one has to introduce the indefinite metric. The commutation relations for 
creation, annihilation operators of the massive state have the wrong sign (— ) 



4 



[2] . The Hilbert space of states is the direct sum of the two subspaces H p and 
H n with positive (H p ) and negative (H n ) square norms. The massless states 
correspond to a positive square norm, and massive states — to a negative 
square norm. The transitions between two subspaces H p and H n break the 
unitarity of the theory. But if the mass m — )■ oo such transitions are forbidden 
and the unitarity is recovered. Thus, a ghost can be removed in the theory 
at large m. This procedure is similar to the Pauli-Villars regularization of 
Feynman diagrams. Therefore, there is physical sense of the Podolsky theory. 
We also can argue (similar to Lee- Wick model [8]) that there is no problem 
with unitarity if the massive photon decays to ordinary fermions through its 
couplings and is not in the spectrum. 

2.2 First-order field equations 

Now we reformulate the third-order field equation (2) in the form of first- 
order relativistic wave equation. Let us consider the system of first-order 
equations 

<9 M W + m$ v = 0, (6) 

dutpn - dp^v + rra/v = 0, (7) 

+ = 0, (8) 

ciVV _ di^v + rmjjpv = 0, (9) 



where 



1 1 

iftfj, = mAp, ip^ = F^, ip p = —d v F vll , ip^ = —d^F^. (10) 



After replacing %p u from Eq.(6) and ip pv from Eq.(9) into Eq.(8), one obtains 
Eq.(2). Eq.(7) is the usual equation for the potentials. Thus, we claim that 
the system of first-order equations (6)- (9) is equivalent to the third-order 
Eq.(2). Let us introduce the 20- dimensional wave function 



( M*) \ 



(A = ti,\fiv],ji,\fiu]), (11) 
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where = Vv( x )> V^(^) = fcjfc) = The function 

represents the direct sum of two four- vectors ipn(x), ipfj,(x), and two 
antisymmetric tensors of the second rank ip^i^x), i/j U v(x)- 

We explore the elements of the entire matrix algebra e A)B [18] . [T9] with 
matrix elements and products 

f^M,N\ r r -M,A^B,N r _M,JV /-irA 

(e J = OmaOnb, e ' e ' = o A b£ , (12) 

where A, 5, M, N = fi, [fiu],Jl, [/xz/], and generalized Kronecker symbols 

= 8 v /3 — Seduce- 

The e M,Ar are 20 x 20-matrices, that consist of zeros and only one element is 
unity where the row M and the column N cross. 

With the help of Eq.(ll),(12) the system of equations (6)-(9) can be 
represented in the form of the first-order equation 

(V' M + e [v ^ v + + e^A m B (x) 

V / AB 

(13) 

+m fi e M.M + e v $ + e u > v + - £ ^ V A * B (x) = 0. 
V 2 2 / ab 

There is a summation over all repeated indices. We define 20-dimensional 
matrices as follows: 

/3 U = pjp + pjp = e v ^ + e^' v , = e"^ + e [u ^, (14) 

p = i^M-M _|_ _|_ _j_ 1 £ M,M_ 

Taking into account Eq.(14),(15), Eq.(13) takes the form of the first-order 
relativistic wave equation: 

(^ + mP)*(i) = 0. (16) 

The presence of the projection operator P in Eq.(16) is connected with the 
fact that there is a massless state in the spectrum [19], [20]. Thus, we re- 
formulated the higher derivative equation (2) in the form of the first-order 
equation (16). The P is the projection operator, P 2 = P [21] and it is not 
the Hermitian matrix P + ^ P. The matrices (3^ and (3^ are Hermitian ma- 
trices and have non-zero components in 10-dimensional subspaces (//, [/if]), 
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(/x, respectively and obey the Petiau-Duffin-Kemmer algebra [22], [23] 

(see also [IB], [IS]): 

PfiPuPa + PotPvPp, = SpwPa + SavP/i- (17) 

Therefore, the matrix 0^ is the direct sum of two 10-dimensional Petiau- 
Duffin-Kemmer matrices. The projection operator P "connects" two 10- 
dimensional subspaces (/x, \jjlu]), and (/J, \jjlu]). Thus, HD Podolsky's electro- 
dynamics equations lead to "doubling" the dimension of the Petiau-Duffin- 
Kemmer algebra representation. 



3 The Lorentz covariance and Hermitianizing 
matrix 

Let us prove the Lorentz covariance of Eq.(16). The Lorentz group transfor- 
mations of coordinates are given by x'^ = L^ v x' v , where the Lorentz matrix 
L = {Lyv} satisfies the equation L m L VOL = 8^ v . The wave function (11), 
under the Lorentz coordinates transformations, becomes 

#V) = T^{x), (18) 

where 20 x 20— matrix T realizes the reducible tensor representation of the 
Lorentz group. The first-order wave equation (16) is transformed into 

(P^ + mP) ¥(x') = {jP^dy + mP) TV(x) = 0, (19) 

where d'^ = L^ u d v . We have the Lorentz covariance of Eq.(16) if equations 

(3^TL^ = T(3 UJ PT = TP (20) 

hold. The infinitesimal Lorentz matrix is given by the 

f J fiv $fj,u ~\iv; ~flf ~v\i: (21) 

where are six parameters defining rotations and boosts. The matrix T 
at the infinitesimal Lorentz transformations reads 

T = 1 + -e^J^, (22) 
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where J^ v are generators of the Lorentz group in 20— dimensional space. 
With the aid of Eq.(21),(22) (using the smallness of parameters e^), we 
obtain from Eq.(20) 

fifiJav Javfifi ^a^il^v ^ufiPay PJav J avP ■ (23) 

The Lorentz group generators in the 20- dimensional representation space are 
given by 

— £ n> v _ e v,i* _|_ £ i x fA,[M _ e M,M (24) 

— 4- ^IVLM _ e [Aw],[A|t] 

and obeys Eq.(23). Thus, we have proved the Lorentz covariance of first-order 
wave equation (16). In Appendix A, we generalize equations considered on 
the case of field equations with the source. It is easy to verify with the help 
of Eq.(12) that the generators (24) obey the usual commutation relations 

[Jfivi Jap] — + Sfi^Jua ~ K^J^a — <^a^/3- (25) 

The Hermitianizing matrix rj should satisfy the relations [24] 

V0m = -P+V + , vPi = Ptv + (m = 1,2,3). (26) 

We find 



^ g.m,m ^A-A _|_ ^[m4],[m4] _^[mn],[mTi] 

_|_^.m,m £^j4 -X- g-[ m 4],["i4] ^[mti],[mti] 



(27) 



2 

The matrix rj is the Hermitian matrix, rj + = rj and commutes with the 
projection operator P: 

v p = p v . (28) 
Consider the "conjugated" wave function 

*(a?) = V + (x)ri = ^, , (29) 

and ty + (x) is the Hermitian conjugated wave function. We took into account 
that for neutral fields, (Vv^^o) are rea l variables. Thus, the relativistically 
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invariant bilinear form is ^/(x)^/(x) = ^/ + (x)i]^/(x). Then, we obtain from 
Eq.(16) the "conjugated" equation 



V(x) 



ix ^ H 



mP + ) 



Formally, one can construct the Lagrangian 
1 



C —2 



0. 



(P^ + mP) m{x) - V(x) - mP + ) *(x)] . (31) 



(30) 



By varying the action S = J d 4 xC, corresponding to the Lagrangian (31), we 
obtain equations of motion (16), (30). One can check using Eq.(26),(29) that 
the Lagrangian C is the real function, C* = C In addition, for neutral fields 
the equality 

^(x)P + ^(x) = V(x)PV(x). (32) 

is valid although P + ^ P. If one wants to consider charged fields (not photon 
fields), then the electric current density is given by 



(33) 



were the complex conjugation * does not act on the metric imaginary unit. 
Using equations of motion (2), one can verify that electric current is conserved 
d^J^x) = 0. The electric current density (33) can be cast into the matrix 
form 

J li (x)=M(x)pW*(x). (34) 

It follows from Eq.(33) that for the neutral (photon) fields the electric current 
density vanishes, J^(x) = 0, as it should be. 



4 The mass and spin projection operators 

Let us consider solutions to Eq.(16) with definite energy and momentum. In 
the momentum space Eq.(16) becomes 

A*(p) = 0, A = ip + mP, p = (3^, (35) 

where p^ is a four- momentum p^ = (p, ipo). Let us consider the massive state, 
p 2 = — m 2 . For this case the 20-dimensional matrix A obeys the equation 
(see (B5) in Appendix B) 

A (A - m) (A - 2m) (A 2 -mA- m 2 ) = 0. (36) 
9 



From Eq.(36), we find the solution to Eq.(35) in the form of the matrix 

II = iV (A - to) (A - 2m) (A 2 - toA - to 2 ) , (37) 

where N is a normalization constant, so that All = 0. This means the the 
every column of the matrix II is the solution to Eq.(35). The requirement 
that the II is the projection operator, II 2 = II, leads to the normalization 
constant N — — 1/ (2m 4 ) [21]. The projection operator (37) extracts solutions 
to Eq.(35) for definite energy and momentum corresponding to the massive 
state. 

With the help of the Lorentz group generators (24), we obtain the spin 
operator (see [2Tj): 

% % 

a P = ~ 771 i e abcPaJbc = — 1 ^abcPaPbPc- (38) 

2 | p | | p | 

The operator (38) obeys the "minimal" matrix equation: 

a p (a p - 1) (<7 P + 1) = 0. (39) 

In accordance with the general method [21], we obtain the projection oper- 
ators extracting spin projections ±1 and 0: 

S(±i) = ^(T P (cr p ± 1) , % = 1 - a 2 , (40) 

satisfying the relations: Sh.^ = £(±1), S^ii^o) — 0, Sf \ = 5( ). 

One may check with the help of Eq.(12) that the operators (40) commute 
with the mass projection operator (37). As a result, from Eq.(37),(40), we 
find projection operators 

A ±1 = ILS (±1) , A = nS (0 ) (41) 

extracting solutions to Eq.(35) for definite energy-momentum, spin projec- 
tions ±1, 0, for states of particles with the mass to. Eq.(41) defines also the 
density matrix for pure spin states. It follows from "minimal" polynomial 
equation (B4) that for the massless state, p 2 = 0, zero eigenvalues of the ma- 
trix A are degenerated, and therefore it is imposable to construct solutions 
to Eq.(35) in the form of the projection operator [21 J . 
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5 Quantum mechanical Hamiltonian 



Now we obtain the quantum mechanical Hamiltonian from equations (6)- 
(9). The Schrodinger form of equations has some attractive features because 
non-dynamical components of the wave function are absent. To find the 
Schrodinger form of Eq.(6)-(9), we exclude the non-dynamical components. 
Eq.(6)-(9) can be cast in the form of two systems 

mip 4m = d 4 ip rn - d m ip 4 , m^4m = d^ m - d m $ 4 , 

(42) 

d 4 ip m4 + d n ip mn = -mif) m , d 4 i(j m4 + d n ip mn = -mijjm, 

m^mn = d m 4> n - d n ip m , mtjj mn = d m tjj n - d n tjj m , m$ 4 = d m $ m4 . (43) 

We can to exclude auxiliary (non-dynamical) components ip mn , ipmn, ip4 from 
Eq.(43). However, the ip 4 can not be excluded from Eq.(42). To introduce 
the evolution of the ip 4 in time, we use the Lorentz condition d m ip m + d 4 ip 4 = 
0. After the replacing the non-dynamical components ip mn , ip mn , i>4 from 
Eq.(43) into Eq.(42), we obtain the equations as follows: 

id t ip m = m^ m4 - d m ip 4 , 

id t 4> 4 = d n ip n , 
idt^m = m^jj m4 - d m ^ 4 , 

l4>n + — (d m d n 1p m - d^lpn) , 

m v ' 
mil> n + — (d m d n i) m - d^tpn) . 



(44) 



id t ip. 



n4 



mi 



id t ip. 



Eq.(44) show that 13-components of the wave function ^(x) possess the 
evolution in time. Therefore, we introduce the 13-component wave function 



( M x ) \ 

1prn 4 (x) 



(45) 



With the help of the elements of the matrix algebra Eq.(12), we rewrite 
Eq.(44) in the Schrodinger form 



(46) 
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where the Hamiltonian is given by 

U = m ( e n ' [ni] + e™' [n4] + 5 [n4] '™ + e [n4] 'A + (e 4 ' m - e mA ) d n 



1 

+— 

m 



£ [m4],n + £ [m4],n _ £ n,[m4]\ _ ( £ [mA],m + £ M], m ^2 



(47) 



From the minimal equation (C6), obtained in Appendix C, we find the pro- 
jection operators extracting states with positive and negative energies for the 
massless states (p 2 = 0) 



o _ (n± |p|) n 2 {U 2 - p 2 - m 2 ) (U 2 - 2p 2 - m 2 ) 



El = ± 



2|p| 3 m 2 (p 2 + m 2 



(48) 



and massive states (p 2 = —m 2 ) 



(H±p )H 2 (U 2 -p 2 ) (U 2 



Pi) 



(49) 



2plm 2 p 2 

Projection operators (48), (49) obey equations as follow: 

(s° ± ) 2 = s° ± , nTP ± = ±p Q TP ± (po = IpI) 

(S ± ) 2 = E ± , HZ± = ± Po Z ± (p = v/|p| 2 + 

Projection operators (48), (49) can be used to construct physical states in 
13-dimensional space of wave functions (45). 



nr 



(50) 



6 Conclusion 

We have formulated Podolsky's generalized electrodynamics equation with 
higher derivatives in the form of the 20-component first-order relativistic 
wave equation. This equation describes vector particles possessing the phys- 
ical massless state and the massive state that is a ghost. To have the con- 
sistent theory the mass of the vector state should be very large. One can 
speculate that the massive vector particles can be described in the gauge- 
invariant manner by this theory. To have the massive state to be the physical 
state, we have to use the reverse-sign in the Lagrangian. Then the Hamil- 
tonian also changes the sign. In this case, however, the massless state be- 
comes the ghost and the question arises: how to get rid of it? Therefore the 
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description of massive particles by Podolsky's generalized electrodynamics is 
questionable. The relativistically invariant bilinear form, and the Lagrangian 
were obtained, and these allow us to use the advantages of the formulation 
of relativistic wave equations. The density matrix obtained can be used 
for quantum electrodynamics calculations in the first-order formalism. It 
should be noted that the Petiau-Duffin-Kemmer form of equations was used 
in quantum chromodynamics i.e. in non-Abelian theory. 

The 13 x 13-matrix Schrodinger form of the equation is derived, and 
the Hamiltonian is obtained. We found projection operators extracting the 
physical eigenvalues of the Hamiltonian. The Schrodinger picture has some 
advantages by considering field interactions. 



Appendix A 



Let us consider the field equation (2) with the source of electromagnetic 
fields — the charge current density: 



(Al) 



We have introduced the current j v with the same dimension as in classical 
electrodynamics. The first-order equations (6), (7), (9) reman the same but 
Eq.(8) is replaced by 



<V/V + ™/v = 3v{x). 
Then Eq.(16) becomes 

(f3 tl d fi + mP)V(x) = P Q j(x) 



where 



3{x) 



( j/x) \ 

jfj,v(%) 



and Pq is the projection operator, Pq 



(A2) 



(A3) 



(A4) 



Pq. The projection oper- 



ator P extracts only the current j M . Therefore, the currents j^x), j^ u (x), 
jixu(x) do not present in the theory and can be put zero. At the Lorentz 
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transformations, j'(x) = Tj(x), and the Lorentz covariance of Eq.(A3) fol- 
lows from Eq.(20),(23) and 

P T = TP , P^Jfiv = J^uPq- (-45) 

The Hermitianizing matrix rj (27) commutes with P , i]P = P rj. Then 
Eq.(30) is replaced by 

- mP + ) = j(x)P , (A6) 
where j(x) = (j^x), -j^(x), j»(x), 

~jfii/( x ))- We obtain the classical limit 
at to — > oo (a — > 0) for Maxwellian electrodynamics from Eq.(Al): 

d,F VfM = j u . (A7) 

Thus, Eq.(A7) is the standard Maxwell equation with the source term. 

Appendix B 

With the help of Eq.(12), we obtain products of matrices entering Eq.(35): 
p 3 = p 2 p i pP + Pp = PpP + p, p 2 P = Pp 2 ) (Bl) 

pPp 2 =p 2 pP, pPpil - P) =p 2 {\ -P), (1 - P)p 2 P = 0. (B2) 
Using Eq.(14), (Bl), (B2), one finds 

A (A — m) = imPpP — p 2 , 



A (A - m) A (A - to) + 2p 2 (A - to) + rap 2 = -ip 4 p - mp 2 p 2 P, (B3) 
A (A - to) A (A - to) - to (A - to) + 2p 2 = imp 2 p - p 2 p 2 - m 2 p 2 (1 - P) . 

From Eq.(Bl)-(B3), we obtain "minimal" polynomials of the matrix A for 
two states: 

A 2 (A - mf = 0, p 2 = 0, (BA) 

A (A -to) (A -2m) (A 2 - mA - m 2 ) = 0, p 2 = -m 2 . (B5) 

It should be noted that zero eigenvalues of the matrix A for the massless 
state are degenerated. 
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Appendix C 



From Eq.(47), we obtain the Hamiltonian in the momentum space 
U = m (e n ' [ni] + e"' [n4] + e [n4] '" + e^A + ip m (e 4 ' m - e mA ) 



1 

+— 

m 



£ [mA],m + £ [m4],m j p 2 _ [ £ [m4],n + £ [m4],n _ £ n,[m4] ( 

Using Eq.(12), one finds 

- p 2 = m 2 f e M-M + gM-M + + e ; 



• (ci) 



% 2 - p 2 - m 2 = m 2 (e 



.[n4],[n4] _ _[n4],[n4] 



(C2) 



(C3) 



Multiplying Eq.(C2) and Eq.(C3), we obtain 

(?/ 2 - p 2 ) {U 2 - p 2 - m 2 ) = (m 2 + p 2 ) p m (e A ^ - e A ^ 
+ (m 2 + p 2 ) p mPn ^M,M] _ £ [n4],M]^ (C4) 
Squiring Eq.(C4), one finds 

(^ 2 -p 2 ) 2 (^ 2 -p 2 -m 2 ) 2 

= p 2 (m 2 + p 2 ) (?/ 2 - p 2 ) (u 2 - p 2 - m 2 ) . (C75) 
From Eq.(C5), we obtain the "minimal" polynomial of the Hamiltonian 

U 2 (n 2 - p 2 ) {U 2 - p 2 - m 2 ) (H 2 - 2p 2 - m 2 ) = 0. (C6) 

Eigenvalues of the Hamiltonian squared read from Eq.(C6): pi = 0, p 2 , = P 2 > 
Po = P 2 + m2 ; = 2p 2 + rn 2 ■ Thus, there are two physical eigenvalues, 
Po — P 2 ; Po — P 2 + m2 > corresponding to massless and massive states of the 
field, and two nonphysical eigenvalues. Eq.(C6) can be used to find projection 
operators extracting physical states in the Schrodinger picture. 
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